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A formula suitable for a quantitative evaluation of the tunneling effect 
in a ferromagnetic particle is derived with the help of the instanton method. 
The tunneling between n-th degenerate states of neighboring wells is domi- 
nated by a periodic pseudoparticle configuration. The low-lying level-splitting 
previously obtained with the LSZ method in field theory in which the tun- 
neling is viewed as the transition of n bosons induced by the usual(vacuum) 
instanton is recovered. The observation made with our new result is that the 
tunneling effect increases at excited states. The results should be useful in 
analyzing results of experimental tests of macroscopic quantum coherence in 
ferromagnetic particles. 
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I. INTRODUCTION 



Macroscopic quantum effects in magnetic systems are of considerable interest both the- 
oretically and experimentally 0. In the context of these investigations the usual termi- 
nology is that macroscopic quantum coherence || (MQC) refers to the resonance between 
neighboring degenerate wells. Some years ago it was reported that MQC was observed for 
antiferromagnetic particles M in resonance experiments. However, the interpretation is con- 
troversial Physically the result of an earlier resonance experiment on ferromagnetic 
particles is not clear since fundamental discrepancies remain between the experimental data 
and theoretical expectations on the basis of magnetic quantum tunneling [|5]||. Apart from 
some other reasons |5|f|, which hinder the acceptance of the observation as definite proof 
of MQC, there is an essential difficulty related to the existing theory of quantum tunneling 
itself in the absence of an external magnetic field. The difficulty was pointed out in ref. @: 
The argument of the WKB exponential of the tunneling for a ferromagnetic particle is 2a/As 
with A = Jp, K\ and K 2 being the hard and medium axis energies. Since s ~ 500 to 5000, 
unless |^ <C 10~ 4 , the tunneling frequency u c is expected to be unobservably small. 

It seems that the unobservably small effect of tunneling at the vacuum level is a common 
phenomenon in various problems, as for instance, in the case of baryon- and lepton-number 
violation at high energy and in the case of pair production of black holes in quantum gravity. 
The one loop correction which results in a prefactor of the WKB leading order exponential 
does not enhance the tunneling significantly in this case. It is a natural and interesting 
question to ask whether the tunneling effect is enhanced by considering tunneling at the 
level of excited states. However, the instanton method is suitable only for the evaluation 
of the tunneling effect at the vacuum level, since the usual(vacuum) instantons satisfy the 
vacuum boundary conditions. Motivated by the study of baryon- and lepton-number viola- 
tion at high energy, recently new types of pseudoparticle configurations were found [f?HTT 



which satisfy periodic boundary conditions and are called periodic instantons |||| or nonva- 
cuum instantons [TOP] • These periodic instantons have, for instance, been used to evaluate 



2 



quantum tunneling at high energy |12|| . There it was confirmed that in the low energy re- 
gion the tunneling effect indeed increases exponentially with energy JI!^| . This finding can 
be expected to have its correspondence in the theoretical analysis of MQC In the present 
paper we therefore adopt the periodic instanton method in order to calculate the tunneling 
amplitude between asymptotically degenerate excited states. 

We derive a compact formula for the level-splitting induced by tunneling which is valid 
for the entire region of energy The results of the application of a method previously 
developed for the calculation of tunneling effects at the level of excited states and based on 
the LSZ (Lehmann, Symanzik and Zimmermann) procedure of field theory |14| are recovered 
in the low energy region jnj. In particular our formula agrees exactly with the level splitting 
of the ground state obtained by means of the usual instanton method [|T6| , |T7| . It is remarkable 



that the tunneling effect enhances significantly if tunneling at the level of an asymptotically 
degenerate excited state is considered. We have shown elsewhere that in certain restricted 
parameter domains the leading contributions of the effect can also be obtained much more 
easily with Schrodinger quantum mechanics |T9| , even in the presence of an applied magnetic 



field 



II. THE EFFECTIVE LAGRANGIAN WITH THE PERIODIC POTENTIAL, 
AND THE ENERGY SPECTRUM FORMULA 

We begin with the following operator Hamiltonian of the ferromagnetic particle which 
has been the starting point of numerous investigations 

H = K 1 s 2 z + K 2 s 2 y (1) 

and describes XOY easy plane anisotropy and an easy axis along the x direction with 



Ki > K 2 > 0. In eq.(|I]) Sj, i = x,y,z, are spin operators obeying the usual commutation 
relation [sj, Sj] = ieijkSk (using natural units throughout, i.e. h — c — 1). Starting from the 
coherent state representation [BTJ of the time evolution operator with Hamiltonian given by 
eq.(|l]) and with the help of the coherent state path-integral we obtain 



( n/ | e - 2 ^ r | ni ) = e^*f-**>K^ f ,t f ;M (2) 

where 

Kfo, t f - U) = J V<pVpe^ C(Mdt (3) 



is the path integral in phase space with canonical variables (f) and p = scos# ||21|| . Also 

C = j>p-H(4>,p) (4) 
is the phase space (or first order) Lagrangian. The Hamiltonian 

has position dependent mass m((p) and potential 

"W ^fl-Wfl - n*) = W' + Q«*+ (6) 

respectively where A = The kets |rij > and |n/ > denote the initial and final spin- 
coherent states and tf — U = IT denotes the difference of final and initial times. Here 
s — s(sin#cos<^,sin#sin0, cos#) is visualized as a classical spin vector with spin number s, 
polar angle 9 and azimuthal angle (p. In the above derivation, the large spin limit s ^> 1 
has been used since giant spins with spin quantum number s ^> 1 are believed to describe 
ferromagnetic grains. A novel feature of the transition amplitude given by eq.(0) is the phase 
factor e~ l ^ f ~^ s which can be put into the Lagrangian, i.e. the expression 0/ — 4>i — It! 
and identified as a Wess-Zumino term |22|. Integrating out the momentum in the path 



integral eq.([|), we obtain the usual Feynman propagator in configuration space, i.e. 

K{4> f ,tfiM= fv<j>e ir H C ^ )dt (7) 



where T>(p is the measure- modified functional differential resulting from the ^-dependent 



M-X 



n - n ^* 



k " 27TZe 



with the second order Lagrangian 

C = imfo) 4? - V{4>) (8) 

which is more convenient for the instanton method used in the following. The potential 
V{4>) is periodic with period n (Fig. 1) and there are two minima in the entire region 2n. 
We may look at this periodic potential as a superlattice with lattice constant 7r and total 
length 2ir, and we can derive the energy spectrum in the tight-binding approximation. The 
translational symmetry is ensured by the possibility of successive 2n extensions. 

It was shown in a previous paper that if e m are the degenerate eigenvalues of the 



system with infinitely high barriers, the energy spectrum is given by 

Em = e m - 2Ae m cos(s + £)vr (9) 

where the expression Ae m is defined by 

Ae m = -[ u* m (<f>, <$> n )Hu m {(f>, $ n+1 )d<p (10) 

which is the usual overlap integral or 2Ae m simply the level shift due to tunneling through 
any one of the barriers. Here u m (4>— $„) denotes the eigenfunction of the harmonic oscillator- 
approximated Hamiltonian Hq in the n-th well, i.e. 

p 2 1 
2m 2 

with m = an d ujq 2 = 4KiK 2 s(s + 1). £ is an integer and here can assume only either of 
the two values "0" and "1". For half integer spin s the spectrum eq. @ is quenched to a 
single degenerate level with degeneracy two. The quenching is seen to be a consequence of 
Kramer's theorem which says that for half integer spin s the degeneracy cannot be removed 



in the presence of the crystal field p3] 
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III. LEVEL SPLITTING OF GROUND STATE DERIVED WITH THE USUAL 

INSTANTON METHOD 



Above we recalled the low-lying energy spectrum of a giant spin particle in the large spin 
limit as obtained in [P3fl . The energy spectrum eq.(|S]) leads to the level splitting given by 
the absolute difference for£ = and 1, namely 



The only parameter left unknown is the overlap integral or the level shift 2Ae m which 
can be evaluated with the help of the instanton method. Instantons in field theory of 
1+0 dimensions are viewed as pseudoparticles with trajectories existing in barriers, and 
are therefore responsible for tunneling. Since instanton trajectories interpolate between 
degenerate vacua (see Fig. 1 (a)) and satisfy vacuum boundary conditions, the instanton 
method is only suitable for the calculation of tunneling splitting 4Aeo between neighboring 
vacua. In the following we first consider tunneling at the vacuum level (i.e. m = 0), which 
leads to the level shift of the ground state energy, i.e. 2Ae of eq.fllOD. Passing to imaginary 
time by Wick rotation r = it, (3 = iT, the amplitude for tunneling from the initial well, 
say that with n = (and fa — 0), to the neighboring well with n — 1 (and <pf — tt) and 
considering large /3, the amplitude for the transition between the corresponding coherent 
states can be shown to be 



AE m = 2Ae m |cos(s + l)n — coss7r| = { 



4Ae, 




for integer spin s 
for half integer spin s 



(11) 



< n(7r)|e- 2/3H |n(0) > = < n(7r)|0,$i >< 0,$ |n(0) > e" 2/3e ° sinh(2/3Ae ) 



7r,P;fa = 0,-P) 



(12) 



where 




is the Euclidean propagator with Euclidean action defined by 




(13) 
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A relation similar to eq.(fL2D applies for a transition between asymptotically degenerate states 
"m" with e Q replaced by e m . Such transitions will be considered in the next section, and the 
splitting will similarly be read off from eq.(^). From now on = denotes the imaginary 
time derivative. 

In the following the Euclidean propagator JCe is evaluated with the instanton method. 
After evaluation we compare the result with eq.(|l4|) to find the level shift 2Aeo- The 
instanton configuration which minimizes the Euclidean action Se is 

<p c = arcsinfcosh 2 uj (t — r ) — A sinh 2 lu (t — t )]~ 1 (14) 

with position r . The instanton trajectory is shown in Fig. 1 (a) with r = 0. The Euclidean 
action evaluated for the instanton trajectory eq. flTjj) , sometimes called the instanton mass, 
is 

/oo . , i _i_ »/Y i i i »/Y 

_ m m^Mldr = y/s(s + l)h. « (s + -) In (15) 



in agreement with refs. [ IB| , r7 | and with refs. [0,0 in the case of small values of A. To our 
understanding it is argued in refs. [|i"6i|T7f that except in e~ Sc one can replace s(s + 1) by 



s 2 , whereas in the exponential factor ys(s + 1) has to be approximated by s + |, even for 
large s. The functional integral K,e can be evaluated with the stationary phase method by 
expanding about the instanton trajectory cj) c such that cj) = <fi c + r], where rj is the small 
fluctuation with boundary conditions r)(/3) = i](—(3) = 0. Up to the one-loop approximation 
we have 

K E = e~ Sc I (16) 



where 



/= / V V e- 5SB 17 

M-0)=o 



is the fluctuation functional integral with the fluctuation action 



SS E = [ r]Mr]dr (18) 
J-p 



where 



M = -^M0 c )|: + Wc) (19) 

with 

V{(j) c ) = \[-m'{<f> c )^ c - iro"(0 c )# + y"(0 c )] (20) 

Here M(4> c ) is the operator of the second variation of the action and m'((f) c ) = dn ^ \et>=<p c - 
As in the usual method of evaluating the fluctuation integral /, we expand the fluctuation 
variable rj in terms of the eigenmodes of M and set rj = H n C n tp n , where ip n denotes the n-th 
eigenfunction of M, and express the result of the integration as an inverse square root of the 
determinant of M. In view of the time translation symmetry of the equation of motion, the 
functional integral K,e is not well defined when expanded about the classical solutions. The 
translational symmetry results in zero eigenmodes of the second variation operator M of the 
action (which in the present case, of course, has only one). This problem can be cured by 
the Faddeev-Popov procedure [[23J] or in a more systematic way with the help of the BRST 
transformation [|25|| . Following the procedure of refs. P^PS] the one instanton contribution 



to the propagator in the one-loop approximation is calculated to be 

/Cg) = 2/5-(l - \Y^s 2 K 2 e- UoP e- Sc (21) 

7T 

To obtain the desired result proportional to sinh(2/?Ae ) in eq. ([12]) , the contributions of the 
infinite number of instanton and anti-instanton pairs to the one instanton contribution have 
to be taken into account (the trajectory of one instanton plus a pair is shown in Fig. 1 (b)). 
Interactions among instantons and anti-instantons are neglected in the dilute instanton-gas 
approximation. The contribution of one instanton plus n such pairs to the propagator is 
obtained with the help of the group property of the Feynman path-integral and is found to 
be 

4 2n+1) = ¥r^S-) hx ^ 22 {T^^ ( 22 ) 

[Zn-\-l)l % [I — A)TV 



Summing over all contributions K,^ n+1 \ the final result of the propagator is found to be 

K E = Ai(-)3e-^° s inh[2/5.2 2 {-^-} iAisf e " & ] (23) 

7T (1 — A)7T 



Comparing with eq.(|T2]) the level shift is seen to be 2 Aeo (note that this is the shift a single 
level) with 

Ae = (24) 

(1 - A)tT 

Expanding Aeo in the region of small values of A, this is 

Ae = 2 2 {^}hMe- 2 ^ (25) 



7T 



the level splitting being 4Ae . Since s is a large number (500 to 5000 as cited in the literature 
@), the level shift is suppressed in leading order by the factor 

e -s c _ e -2sV\ (26) 

For Aeo to be observable, plausible values of A are of order 10~ 5 . Such small values of A may 
not be realizable in nature. The quantum correction proportional to ei lns does not increase 
the tunneling probability significantly. 

It may be interesting to compare our level splitting of the ground state AE with that 
in refs. [^,0. The Hamiltonian of eq.(l) in refs. and ]17| equals that of our model if 



one sets the external field equal to zero (h = 0) and makes the replacements A = K 2 , B = 
Ki — K 2 . By simple algebra one can find that the level splitting eq.(16) in ref. [|17| (or eq.(9a) 
in ref. |Tj|) coincides exactly with our AE = 4 Aeo which also simplifies to the results of refs. 
19| , |20|p3| in the appropriate limits. We can now proceed to consider nonvacuum or periodic 



instantons which are the pseudoparticles interpolating between asymptotically degenerate 
excited states, and which reproduce the above result for vacuum instantons in the vacuum 
limit. 
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IV. QUANTUM TUNNELING FOR EXCITED STATES AND THE 
GENERALIZED FORMULA FOR THE LEVEL SPLITTING 



The periodic instanton method @-12] has become a powerful tool for the evaluation of 
quantum tunneling amplitudes over the entire region of energy. The model at hand can be 
looked at as one for tunneling at the level of excited states of a sine-Gordon-type potential 
with a position-dependent mass,which has not been reported previously in the literature. 
The periodic instanton configuration which minimizes the Euclidean action eq. ( |T3"D is seen 
to satisfy the equation of motion 



m(4> p ) ( c% 



dr 



V{<j> p ) = -E c 



(27) 



where E c i > 0, which is a constant of integration, may be viewed as the classical energy of 
the pseudoparticle configuration. Through the usual procedure of derivation of the periodic 



instanton solution, after a laborious but straightforward integration of eq. (27Q, we obtain 



the periodic instanton configuration [27 



arcsm 



1 — k 2 sn 2 (u>r\k) 



1 — \k 2 sn 2 (ur\k) 
where sn {ur\k) denotes a Jacobian elliptic function of modulus k, 



(28) 



k 



\ 



n 2 



n 2 



A 



(29) 



with 



K 2 s(s + 1) 



U = U Q i 1 



(30) 



E d U V K 2 s(s + 1) 

One can check with sn(w|l) = tanhw that for k 2 = 1 (corresponding to E c i — > 0) the 
configuration (p8|) reduces to the instanton of eq. (|H|) . The Jacobian elliptic function sn has 
real periods An)C(k), n being an integer and JC(k) the quarter period given by the usual 
elliptic integral of the first kind. The parameter values ( |2"9"D , ( |3U| ) ensure the periodicity of 
(p p = at r = —2/5, +2(3 with a crossover from negative to positive values at r = and 
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= ^mod2n. Thus the one-way transition from a turning point a\ to the other turning 
point a2 is mediated by the instanton-like part or one half of the periodic instanton extending 
from r = —f3 to r = +/3 (the periodic instanton itself returning to its original position after 
time 4/3). The trajectory of the periodic instanton eq.(^) is shown schematically in Fig. 
1(c) where the trajectory is shifted by an amount 2tx. The instanton-like part starts at 



time —(3 from turning point <fi = a\ = arcsin y K2S £j +1 ) and reaches the other turning point 



= a,2 = 7r — arcsin J K2S £j +1 j at time (3. The Euclidean action of the periodic instanton 
configuration eq.(p8|) over the domain from r = — (3 to r = +/3 can be found to be (cf. refs. 
10) 



S p 



f [m(0 p )0j + E d ] dr = W + 2E d p (31) 



where 



W 
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JC(k)-(l-\k 2 )Tl(\k 2 ,k) (32) 



and Il(A/c 2 , k) is the complete elliptic integral of the third kind. 

It is now necessary to calculate the amplitude for the transition mediated by one 
pseudoparticle configuration - in the present case the instanton-part of the periodic instan- 
ton - and then to sum over this contribution with other contributions A^ 2n+1 ^ obtained by 
adding an arbitrary number (n = 1,2,3, ...) of noninteracting complete periodic instantons 
(corresponding to addition of instanton-anti-instanton pairs in the dilute gas approximation 
of the instanton procedure) so that in the path integral representation of the amplitude all 
possible paths between —f3 and +(3 are taken into account. As is well known, this results 
effectively in an exponentiation of the first contribution (see e.g. ref. [|ll|]). A calculation 
very similar to this can be found in reference |I(J (cf. there eqs. (4.19), (4.24), (4.25)). 
We can therefore simply transcribe the result into the present context. The total amplitude 
corresponding to eq.(^) is then 

2(3u 



A=J2 A(2m+l) = e~ 2E ^smh 



m=0 



4K,(k' 



_ e -w 
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Comparison with eq. ( 12|) as noted earlier then yields the compact formula [[8] 10 



-w 



AlC(k') 



(33) 



k' being the complementary modulus of k, i.e. k' = yl — k 2 . We emphasize that the above 
formula eq. ( j33l) obtained with the periodic instanton is valid for the entire region of energy 
E d . The validity of the level splitting of excited states considered previously in the literature 
28fl was restricted to the quasi uniaxial limit, 1 — A 1. Our approach, however, is valid 



over the whole range of A. 

We now consider the low energy limit where E c i is much less than the sphaleron energy 
(or barrier height K 2 s(s + 1) of the potential), i.e. k — > l,k' — > 0. Expanding W as 
power series of k' up to quadratic order as in ref. |9j and making use of the oscillator 
approximation of the periodic potential around one of the minima with the quantization 
replacement E d — > e m = (m + |) a> (as in refs. P,|T0|j), we then have 



Vx ( i , , 

+ \ m + - ) In 



1 - A / 1 
— =r- [m H — 

8^ V 2 



m H — 
2 



The level shift 2Ae m is finally given with 

Ae = J_ 2 3m + 2^ (m+1 > ] S m ^ \ K,K 2 



exp 



, i i + Vx 
- (s + 2 )ln TT7x 



(34) 



(35) 



mi (1-A) m Itt(I-A) 

which reduces to the ground state expression eq.(p4|) for m — 0. Eq. (|35D provides the 
linkage from the splitting of excited states (which up to now has been considered only by a 
perturbative method f26|) to that of the ground state. 

It has been pointed out that in the low energy region the level shift 2 Ae m may be obtained 
by an alternative but equally good approach called the LSZ method in which only the vacuum 
instanton (eq. (|14"D for the problem at hand) plays a role [|l^,0,|l^]. The comparison with 
the LSZ method is explained in the appendix. The level shift formula eq. ( |35| ) is suitable for 
a quantitative evaluation of the tunneling amplitude. We therefore borrow from ref. [H] some 
data for a ferromagnetic particle in order to demonstrate the difference in A with respect 
to the result of tunneling at the ground state (but we by no means imply - for the reasons 
explained in the introduction - that the transition is a result of quantum tunneling). 
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Table 1 



Comparison of A— values computed for various excited states 



m 


A for s = 3 x 10 3 


A for s = 2 x 10 3 


A for s = 10 3 


2 


0.17350 x 10" 4 


0.37641 x 10" 4 


1.41226 x 10" 4 


4 


0.24543 x lO" 4 


0.53458 x 10~ 4 


2.01990 x 10~ 4 


6 


0.32066 x 10~ 4 


0.70031 x 10" 4 


2.65869 x 10~ 4 


8 


0.39985 x 10- 4 


0.87500 x 10" 4 


3.33361 x 10- 4 


10 


0.48316 x lO" 4 


1.05899 x 10~ 4 


4.04587 x 10~ 4 


12 


0.57064 x 10" 4 


1.25235 x 10" 4 


4.79555 x 10" 4 


20 


0.98155 x 10~ 4 


3.15062 x 10" 4 


20.29119 x 10- 4 



The theoretical level shift calculated with the simple model of the WKB exponential is 
negligibly small. If the observed frequency reported in ref. || were considered to be the 
result of tunneling from the ground state, A would be 1.01 x 10~ 5 . The tunneling effect 
enhances significantly at excited states. Table 1 shows how A increases with m (the index 
of excited states) for tunneling frequency u c = 500Hz which is the level splitting 4Ae m and 
s = 3 x 10 3 , 2 x 10 3 and 10 3 respectively. For s = 10 3 to attain the frequency uo c = 500Hz 
the value of A need only be 2.03 x 10~ 3 , which may be a physically acceptable quantity for 
a small ferromagnetic particle. 

V. CONCLUSION 

We have shown that the periodic instanton method as well as the LSZ method are useful 
for the calculation of tunneling effects at excited states of a spin system at low temperature, 
the former being valid over the entire region of energy, the latter at low energies. The 
results are suitable for a quantitative evaluation of the tunneling effect. With the simple 
model of eq.([I]) for a ferromagnetic particle we conclude that the tunneling effect at the level 
of excited states increases. 
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Appendix: Comparison with the LSZ result 

From the viewpoint of the LSZ method the transition amplitude between m-th degenerate 
eigenstates in any two neighboring wells (say those for n = 0, 1) is viewed as the transition 
of m bosons induced by the instanton of eq. (|14|) and is related with the level shift 2Ae m by 

A™ = (m, fcile-^lm, %) = e~ 2l3tm sinh2/3Ae m (A.l) 



Following ref. |L3[ the transition amplitude as well as the S-matrix can be related to the 



Green's function through the procedure known as the LSZ reduction technique fll4 |. To this 



end we construct the interacting Euclidean fields <f>± in the classically forbidden region with 
the help of the instanton configuration of eq . (|14"D . Thus we define 



cj) + :=TT- <f> c , </>_ := (j) c (A. 2) 

such that the interaction fields vanish in their respective asymptotic regions, i.e. 

lim 0+ = 0, lim 0_ = (A.3) 



T — >00 T — > — OO 



The subscripts "-" and "+" here denote the wells with minima at $o — and $i = n 
respectively. The Euclidean creation and annihilation operators a± and a± which create and 
annihilate an effective boson in wells "+" and "-" respectively are related to the interaction 
field operators (f)± by 



V UJq OT 



- - -f^°* I few < A - 4 » 
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where m = an d e -g- 



-y/2m uj (l — A) 2 



The transition amplitude of eq. ( |A. 1| ) can now be written 



Am nm — 2/3moJo 



with S-matrix element 



S%= lim — (0|a + (r/)...a + (r/)aL(rl)...aLK)|0) 

■" r*-f-oo m! 



A At 



(A.5) 



(A.6) 



t/- 



The S-matrix element can be evaluated in terms of the Green's function G which arises in 
its evaluation. Thus 



nm 



lim 

rf — ''Oo 



n 

i=i 



UJ 



drf 



u dr\ 



1 m 

in 



e -o(r/-r/) 

^0 



d 2 G fdG dG ^ 

+ ^0 —f ~ -7T1 ~ UqG 



drf drf 



where the 2m-point Green's function is defined as usual, i.e. 



(A.7) 



MO|0> 



(A.8) 



Following the procedure in refs. [12,13,15] the S-matrix element is thus 



S 



771 



where 



A 



2/32 2 



2m c 

LU 



drf dr\ 



4° 

A f,i 



(A.9) 



(A.10) 



:i-a)tt 

is the amplitude for the transition between degenerate ground states. Using the instanton 
solution of eq. fll4l) and taking the large time limit r\ — > — 00,7/ — > 00 after performing 
the imaginary time derivatives, we finally obtain the transition amplitude by observing that 
each pair of vertices in eq. (|A.9|) contributes a factor — 4wq. Then 



A 



1 ( 2 3 m U \ Wonl 2f3 A 



(A.ll) 



Comparing eq.( [A.ll| ) with eq. ( |A.l ) we recover the level shift 2Ae m of eq.(^) for the low- 



lying levels. Following refs. 0)0] we nere take Uq = 2\JK\K2S for reasons of consistency. 
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Figure Captions 

Fig. 1 

The periodic potential and the instanton trajectories: 

(a) For one instanton (i.e. vacuum instanton), 

(b) For one instanton plus one instanton-anti-instanton pair, 

(c) One half of the periodic instanton. 



18 




Fig. 1 



